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Abstract 

We investigate the spin correlations in top quark pair production near thresh- 
old at the e~e + linear collider. Comparing with the results above the threshold 
region, we find that near the threshold region the off-diagonal basis, the optimized 
decomposition of the top quark spins above the threshold region, does not exist, 
and the beamline basis is the optimal basis, in which there are the dominant spin 
components: the up-down (UD) component for e^e + scattering and the down-up 
(DU) component for e^e + scattering can make up more than 50% of the total cross 
section, respectively. 
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The top quark is the heaviest known particle in the standard model (SM) and rapidly 
decays before it hadronizes [1], and the spin information of the top quark is preserved 
from production to decay. Thus we can expect the spin orientation of the top quark to be 
observable experimentally. The spin correlations for the top quark pair production above 
the threshold at e + e~ colliders have been extensively discussed [2]. Parke and Shadmi 
[3] proposed the generic spin basis and found that the "off-diagonal" basis, a special case 
of the generic spin basis, is a more optimized decomposition of the top quark spins for 
e + e~ colliders. As shown in Refs. [4, 5], the off-diagonal basis is indeed the optimal spin 
basis even after the inclusion of 0(a a ) QCD corrections: at ^Js = 400 GeV using the 
off-diagonal basis the dominant spin components in both e~[ and e\ scattering make up 
more than 99% of the total cross section at both tree and one-loop level, but such fraction 
is only ~ 53% in the helicity basis. 

Up to now, all studies of spin correlations of the top quark pair production at the 
e~e + collider are limited only to the process of ti production above the threshold region. 
In the future e~e + linear collider, threshold production of the top quark pair will allow to 
study their properties with extremely high precision. Because of large top quark mass and 
decay width, the bound-state resonances lose their separate identify and smear together 
into a broad threshold enhancement [6], and as a result, the nonpretubative QCD effects 
induced by the gluon condensate are small, allowing us to calculate the cross section with 
high accuracy by using perturbative QCD even in the threshold region. And it is interest- 
ing to investigate the spin correlations of the top quark pair production near threshold. 
Some methods used to deal with the behavior of top quark pair production near threshold 
have been established. The Green function technique was demonstrated suitable to cal- 
culate the total cross section and to predict the top quark momentum distribution, and 
independent approaches were developed for solving the Schrodinger equation in position 
space and the Lippmann-Schwinger equation in momentum space [6, 7, 8, 9]. 

In this letter, we study the spin correlations in the top quark pair production near 
threshold. In order to include effects of the quark-antiquark potential and the decay width 
in our calculation, we use two ingredients: the first is the vertex function, which represents 
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the QCD binding effects and the anomalous interactions, and can be obtained through 
solving Lippmann-Schwinger equation [10]; the second is the two unstable-particle phase 
space in the nonrelativistic limit [7], which, combined with the vertex function, describes 
the top quark momentum distribution. As a cross check of our calculation, after using 
above ingredients, also we obtain the total cross section and the momentum distribution 
which are consistent with ones given in Ref. [8]. 
In the SM, we consider the process 

e~e + -> it (1) 

at the LC with y/s ~ 2m t . The tree level Ve~e + vertex can be written as 

Kee = nKlP- + K^P + ), (2) 

where P± — (1 ± 7s)/2, and the SM values for these coupling factors are K\ R = — e for 
V — 7, K[ = e(2sin 2 9 W — l)/2 sin9 w cos9 w and = e sin. 0w/ cos9\y for V — Z, and 
the 9 W is the Weinberg angle. 

The Vtt vertex function T vt i (V = 7, Z) can be generally written as 

Ktt = nAyP- + B v P + ) +1 5 ^C v , (3) 

Tfli 

where is the momentum of the outgoing top quark, and Ay, By, Cy are the form factors. 
As shown in Fig.l, the vertex function T v (Ty t i) satisfies following integral equation [7]: 

Ty = X v + j 0- 4 (-4nC F a s )D flu (p-k)YS F (k + ^Ty(k,q)S F (k-^ 1 \ (4) 

where X 1 = 7^ and X z = 7^75, S F is the top quark propagator, and D^ v is the gluon 
propagator. In the nonrelativistic limit, the propagators are replaced by 

-A-KCpOisD^ip) — > iV(p)5^ 5 u0 , (5) 

l+7» _ t± 

S F (k + ^-) — > 1 2 ^ , (6) 

1-7° _ kj_ 

S F (k-h — > % 2 ^ , (7) 



3 



where E^ R = yfs — 2m t . With above approximations, when the QCD binding effects and 
the CP-violating anomalous couplings are included [10], the form factors (A v , B v , Cy) in 
the vertex function T v are given by 

A i = Ul- 2 -^)G(\p t \) V , (8) 

6 71 

((I-^sin 2 M(l-^W t |V 



sin 6 W cos 6V 4 3 ' n 

+ 1 j(l-^)F(\p t \) lf ), (9) 

4 7T 

= 2 -e(l- 2 -^)G(\p t \)<p, (10) 

6 71 

Bz = , . / , ((l-fsin 2 M(l-— MNV 
4 sin 6*14/ cos 6*14/3 7r 

-(1 - ^)F(|p t |V), (11) 
C ex 2 

C 7 = -ied tl (l--^)F(\p t \)<p + -ied tg D(\p t \)<p, (12) 

7T O 

= - , e — a H^(l-^)F(bl|V + ^(^-|sin 2 M^(bllV),(13) 

Sill (7(4/ COS 7T 4 O 
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where d tl , d tz and d tfl are the anomalous couplings of the top quark, and <p 
(E NR + iT t ). The Green-functions (G, F and D) in Eqs.(8)-(13) are the solutions of the 
Lippmann-Schwinger equations [10]: 

(M! _ ( Enr + iF t ))G(E NR , \p\) + / -0^[V(\P- k\)G(E NR , \k\)} = 1, (14) 
(^ - (E NR + ir t ))fiF(ENR, \P\) + / ^r,[V{\p - %\)&F(E NR , \k\)\ = p\ (15) 
- (E NR + iT t ))?D(E NR , \p\) + / j^\V(W- k\)kW(E NR , \k\)} 

which can be derived from Eq.(4), and the QCD potential V(\p — k\) is given by [11] 

V(\p- k\) = V c (\p- k\) + V BF (\p- k\) + V NA {\p- k\) (17) 



with 

vc{\p-k\) = _ 4 ^(ir gn {1 + ( «^z fc i 2 ) )ai + { ^(\p- m f h (18) 

V ' 17 H2 L V 47T 7 V 47T 
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V BF (\p-k\) 



AnC F a s (\p-k\ 2 )[ 



(p x k) 2 1 



1 3i(p*x + 




m 2 |p-£;| 4 4m 2 
((S t + g) ■ (p - k)f 
\p — k\ 2 



2m 2 \p — k\ 2 



(19) 



(20) 



Po ln(|p - £| 2 /A 2 ) /9o ln(|p - £| 2 /A 2 ) 



(21) 



where St and St are the top and antitop spin operators, and 



ai = 43/9, a 2 = 155.842, A = 226 MeV, 



Po = 23/3, ft 



58/3, /z = 20GeV. 



(22) 



Using above vertex functions, we calculate the spin correlations. In the tt center of 
mass frame (CMS), the scattering plane is defined to be the X-Z plane where the electron 
is moving along the +Z direction and 9 t is defined as the scattering angle of the top quark, 
and we also set fa — 0. The Born helicity amplitudes for the process (1) are obtained by 
summing the contributions from both the Z and 7: 



where s = AE 2 is the total energy in CMS, E e is the energy of the electron, and R(s) = 



In the generic spin basis [3] the top quark (anti-top quark) spin states are defined in 
the top quark (anti-top quark) rest-frame, where one decomposes the top (anti-top) spin 
along the direction s t (%), which makes an angle £ with the anti-top (top) momentum 
in the clockwise direction. Thus, the state (ij^) refers to a top with spin in the +s t 
(—St) direction in the top rest-frame, and an anti-top with spin +% (— % ) in the anti-top 
rest-frame. 

In the generic spin basis, the amplitudes M(h e -, h e +, s t , Sf) for the process e~e + — > tt 
can be generally written as 



M(h e -,h e+ , h t , ht) = M(h e -,h e+ , h t , /i f ) 7 + M(h e -,h e+ ,h t , h) z R(s), 



(23) 



s/(s-M 2 )\ 



M(- + t T t T oitfa) = ±2E e [m t (A L + B L ) sin 9 cos £ - (\p t \(A L - B L ) 



At NLC with y/s ~ 2m t , the imaginary part of the Z propagator can be neglected safely. 
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+ cos9E e (A L + B L ))sm£±2E e C L —sm9], (24) 

Tflt 

M(- + t ] t i or tfr) = 2E e [m t (A L + B L ) sin # sin £ ± {E e {A L + B L ) + \p t \(A L - B L ) cos9) 

+ cos£(\p t \(A L - B L ) + E e (A L + B L ) cosfl)], (25) 

M(+ -i T i T ori^) = ±2S c ["it(>lfl + -Bfl)sin0cosf 

+((^-^^1 -cos0£ e (A R + 5 R ))sin£ T 2£ e C^sin0], (26) 

M(+ - t T t^ or i J T ) = 2£ e [m t (A R + sin 9 sin f =p (E e (A R + B R ) + - A R ) costf) 

+(|jT t |(B fl - A R ) + £ e (A fl + costf) cos^], (27) 

where A LR , B L Rl Cl,r are the form factors and defined as 

Al,r = -(Kl R A, + Kl R A z R(s)), (28) 

Bl,r = -(Kl R B, + Kl R B z R(s)), (29) 

Cl,r = -(Kl R C 1 + Kl R C z R(s)). (30) 

Because the Vtt vertex (V=7,Z) has complex structures near threshold, we can not find 
such a spin angle £ that makes Eq.(24) equal to zero, and then there is not the off- 
diagonal basis, in contrast to the case of the above threshold. The amplitudes in the 
helicity basis and the beamline basis can be obtained by setting cos£ = ±1 and cos£ = 
(cosfl + g)/(l + ^cos#) in Eqs.(24)-(27), respectively. 

In the nonrelativistic limit, the differential cross sections of two unstable particle 
production [7] in the generic spin basis can be expressed as 

da(h e -,h e +,s t ,st) ?t f \Pt\ 2 nj/h , . mi2ji-i /qi\ 

d^Te = 8^m!J (E NR -^y + r? mhe -> K+ > St > Si)l dU (31) 

In the numerical calculation, we use the following parameters as standard input [13]: 

a s (M z ) = 0.117, m z = 91.188 GeV, m t = 174 GeV, 

F t = 1.43 GeV, sin 2 9 W = 0.2311. (32) 

We define the fractions of the total cross sections for the different spin components in 
the beamline spin basis as following: 

R(e L ,R, e + , s t , St) = — — -=r . (33) 

VtotaAe L ,R e+ -»■ U ) 
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With E NR = 5 GeV, we have 



R(e L e + 


- hh 


or e R e + 


- hh) 


~ 50%, 


(34) 


R(e~ L e+ 


- hh 


or e R e + 


- hh) 


~ 24%, 


(35) 


R(e~ L e+ 


— ► t^t^ 


or e R 'e + 


- hh) 


~ 13%, 


(36) 


R(e- L e+ 


- hh 


or e R 'e + 




~ 13%. 


(37) 



These results show that in the tt threshold region all spin components can not be ne- 
glected. In Fig. 2 we show the differential cross sections for the precess eJ j R e + — > tt in the 
beamline basis with E^r = 5 GeV. One can see that there is a dominant spin component 
when scattering angle 9 ranges between n/3 and 27r/3. More precisely, according to the 
definition of R in Eq.(33), we integreted the 9 from tt/3 to 27r/3, instead of to tt, and 
then have R(eJ j e + — > t^h 01 e R e+ ~^ ^lh) — 79%. But, as shown in Fig. 3, in the helicity 
basis there are not such dominant spin components. 

Moreover, in our calculation, we considered the Higgs potential effect on the vertex 
functions. Our numerical results show that such effect is very small and can be neglected. 
The anomalous couplings (d tl , d tz , d tg ) are one of several sources to provide CP-violation 
[12], and the numerical results show that these effects on the spin correlations in the tt 
threshold region are very small, too. For example, with taking d tJ = d tz = d tg = 1CT 3 [10], 
the corresponding changes of the spin correlations are smaller than 0.1%. 

To summarize, we have calculated the spin correlations in the top quark pair produc- 
tion near threshold at the e~e + Linear Collider in the SM. We start from the general form 
of the Vti vertex (V=7,Z) near threshold, derive out the amplitudes in the generic spin 
basis, and give the differential cross sections in the NNLO QCD potential. Comparing 
with the previous results above the threshold region in Refs. [4, 5], we find: 

(a) The most important difference between the two regions is that in the 
above threshold region we can find the off-diagonal basis in which only one 
spin component is appreciably non-zero, but in the threshold region the off- 
diagonal basis does not exist. 
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(b) Near threshold the beamline basis is the optimal basis, in which there 
are the dominant spin components: the up-down (UD) component for e~[e + 
scattering and the down-up (DU) component for e^e + scattering can make up 
more than 50% of the total cross section, respectively. 

(c) The observables of the spin correlations near threshold are less advan- 
tageous than ones of above the threshold region. Nevertheless, because of the 
extremely high measurement precision of the top quark pair threshold pro- 
duction, it is still valuable to study the spin correlations in the top quark pair 
production near threshold. 
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Figure 1: Lippmann-Schwinger equation in diagrammatical form. 
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Figure 2: The differential cross sections in the beamline basis for the e LR e + — > tt processes 
with the ift-|-(UU), £j^(DD), £ptj_(UD) and fjij(DU) productions, assuming _Ejv.r = v^ - 
2m t = 5 GeV. 
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Figure 3: The differential cross sections in the helicity basis for the e L e + — > ti^t and 
e^e + — > tL,Rt processes at E^r = s/s — 2m t = — 1 GeV. 
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